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1. Introductory remarks 

Several interesting lattice field theories have a so-called complex action problem, i.e., in some 
parameter range the action S has an imaginary part and the Boltzmann factor exp(— S) cannot be 
used as a weight factor in a Monte Carlo simulation. Examples are lattice QCD at finite density 
or lattice QCD with a theta term. In a quite remarkable series of developments for several simpler 
systems the complex action problem was overcome in recent years by mapping them onto new 
degrees of freedom (dual representation), where the partition sum is a sum over real and positive 
contributions. The dual variables so far mainly consisted of loops that are subject to constraints. 
Such systems can be updated efficiently with algorithms based on the Prokof'ev Svistunov worm 
algorithm [1]. 

In a recent paper [2] we were able to solve the complex action problem for a new class of 
systems: Abelian gauge-Higgs models. In this case the dual degrees of freedom are surfaces for 
the gauge fields which are either open or bounded by flux lines representing the matter fields. In 
this contribution we report on these recent results and present ideas for further development of the 
concepts for Monte Carlo simulations with dual representations. 



2. Abelian gauge-Higgs models as theories of surfaces 

We here consider two abelian gauge-Higgs models, one based on the gauge group Z3, the other one 
on U(l) with two flavors. The conventional action of the Z3 gauge-Higgs model is given by 



z x v<p X,V 



(2.1) 



The first part is the Wilson plaquette action where U XiVp = U XiV U x+ y p U* + ~ y U^ p * s tne plaquette for 
the link variables U X)V E Z3 = {e' 27r / 3 , <?-' 27r / 3 5 1}. The second part consists of hopping terms for the 
Higgs field variables (f) x E Z3 made gauge invariant with link variables. We work on lattices of size 
TV 3 x N t and use periodic boundary conditions for all directions. j8 is the inverse gauge coupling 
and the factor T] controls the coupling between the Higgs- and the gauge fields. Essentially it plays 
the role of an inverse Higgs mass: If the Higgs field is infinitely heavy (t] = 0) it decouples from the 
dynamics of the system. The hopping terms in the 4-direction (= temporal direction) are coupled to 
the chemical potential fl. In the conventional representation the model then has a complex action 
problem at \i > 0. The partition sum is obtained by summing the Boltzmann factor over all possible 
field configurations Z = exp(— S). 

The complex action problem can be overcome completely [2] by mapping the partition sum 
to a sum over new degrees of freedom: Plaquette occupation numbers p x y P E { — 1,0, 1} which are 
attached to the plaquettes of the lattice and represent the gauge degrees of freedom (see also [3] 
for related results in pure abelian gauge theory), and fluxes k xy E {—1,0, 1} assigned to the links 
of the lattice which represent the matter field. The partition sum can be rewritten exactly into the 
form 

Z = C £ £ V P \p,k] tf F [k] W P \p] W F [k] . (2.2) 

{p}« 
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Here C is a constant depending on the parameters j8 and /I, and the sum runs over all configuration 
of the p and k variables. Each configuration is weighted with real and positive factors 



WW = n B i 



Px,vp 



x,v<p 



(2.3) 



where 



B P 
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e 2n + 2 e -n 



eP + 2e~P/ 2 
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(2.4) 
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-1,0,+1 



In addition the configurations are subject to the constraints 



(2.5) 



p:v<p 



p:p<v 



In the two constraints r(n) denotes the triality function which equals 1 if n is a multiple of 3 and 
vanishes otherwise. The first constraint is a product over all sites and at each site forces the total 
matter flux to vanish modulo 3. The second constraint is a product over all links and forces the 
combined plaquette and matter flux at each link to vanish modulo 3. The corresponding admissible 
configurations have the interpretation of closed surfaces of plaquette variabes or surfaces bounded 
by matter flux, and in Fig. 2 we show simple examples of admissible configurations. 

The other gauge-Higgs system we consider is based on the gauge group U(l). In the conven- 
tional representation the fields are the link variables U XjV E U(l) and the two flavors (f = 1,2) of 

(f) 

complex valued scalar fields X V ' with opposite charges. The action reads, 

2 



s = ~P 1 1 Re u x ,v t4+v, P E£+p,v u% +11 

x v<p f=l X 



^)|#| 2 +A(/)|# ) | 4 



(2.6) 



X,V 





Figure 1: Examples of admissible configurations in the dual representation of the Z3 model. 
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Again we use the Wilson-plaquette action for the gauge fields. For the Higgs fields the parameters 
are quartic couplings and = 8 + , with the bare masses. The hopping-terms in 
4-direction are coupled to chemical potentials ji^ and again we have a complex action problem 
when /xW ^0. 

Similar to the Z3 -model it is possible to rewrite the partition sum exactly in terms of a dual 
representation with plaquette variables p for the gauge fields and two sets of flux variables and 
for the two flavors (f = 1 , 2) of matter fields, 

Z = £ W P \p\wi l \k^\l^W^ , (2.7) 

yfe (1) , yfc (2) e (-00,00), /W, /P) e [0,00), p e (-00,00) . 

The additional sets of flux variables appear because for the U(l) case the matter fields also 

have a radial degree of freedom which is absent in the Z3 case. However, it is important to note 
that the link variables for the radial degrees of freedom do not enter in the constraints. As for 
the Z3 case it turns out that in the dual form of the U(l) model all weig hts Wp^ l \w^ 2) (which 
we do not display explicitly here) are real and positive and the complex action problem is solved. 
The constraints for the flux variables are the same as in Eq. (2.5) for the Z3 case, with the 
only difference that the triality function T(n) is replaced by a Kronecker delta S n $ which here we 
denote as 8(n). The constraint ffp that links the plaquette and flux variables reads 

Vp\pM 1) ,*P ) ]=n s ( I [P*,yp-P*-P,vp\- L [Px, P v-p x -p, pv ] + ki]l-k^), (2.8) 

x,v \p.v<p p:p<v J 

where the different signs for and reflect the opposite charges of the two flavors. 

In this study we use thermodynamical observables that can be obtained as derivatives of InZ 
with respect to the various couplings, such as the plaquette ([/), the particle number density n, and 
the corresponding susceptibilities, 

1 d InZ 1 d 2 lnZ _ 1 dlnZ 1 d 2 InZ 

The derivatives can easily be worked out also in the dual representations and the observables are 
then obtained as combinations of moments of the dual variables. 
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Figure 2: Cube updates (two lhs. plots) and plaquette updates (rhs.) used in our local Metropolis algorithm. 
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Figure 3: Examples of segments used in the surface worm algorithm. The plaquette occupation numbers are 
changed as indicated by the signs. The links marked with full (dashed) lines are changed by +1 (—1). The 
empty link shows where the segment is attached and the dotted link is the new position of the defect where 
the constraints are violated. 



3. Update strategies for the dual variables 

In order to do simulations in the dual representation we need a Monte-Carlo algorithm which gen- 
erates admissible configurations of the dual variables, i.e., configurations that obey all constraints. 
The simplest approach is a local Metropolis update where plaquette and flux variables on cubes and 
plaquettes are changed according to a simple pattern, such that the constraints remain intact. We 
use two updates which we refer to as cube and plaquette updates. The cube update only changes 
the plaquette variables p on the faces of a 3 -cube embedded in four dimension according to one of 
the two patterns in the two lhs. plots of Fig. 2, where for the case of Z3 all additions are understood 
modulo 3. It is easy to see that such a change leaves the constraints ffp intact, and since no flux 
variables are changed also the constraints A full sweep of cube updates consists of offering one 
of the two cube updates of Fig. 2 with equal probability to all 3-cubes of the lattice and accepting 
the changes with the usual Metropolis probability. To update the flux variables, in the plaquette 
update we offer a joint update of one plaquette variable and the fluxes on the links of that plaquette 
according to the possibilities shown in the two rhs. plots of Fig. 2 (again addition is modulo 3 for 
Z3). Also here one easily shows that the constraints ffp and ffp remain intact, and a full sweep 
of plaquette updates consists of offering the change to all plaquettes of the lattice. For the case 
of the 2-flavor U(l) model there exist two different plaquette sweeps for the two flavors, and the 
unconstrained variables are updated with conventional Monte Carlo sweeps. 

A second, more efficient strategy is a generalization of the worm algorithm [1] to surfaces 
with boundaries [4]. After changing the flux at a randomly chosen link, which creates a defect 
in the constraints, one adds so-called segments which consist of a plaquette and two of its links. 
The corresponding plaquette number and the fluxes on the two links are changed (see Fig. 3 for 
examples of segments). Attaching such a segment to the defect propagates the defect on the lattice 
until the worm decides to heal the violated constraints at the defect by adding the missing flux. 
It may be shown [4], that the generalized surface worm algorithm clearly outperforms the local 
updates for a wide range of parameters. 



4. Selected results 

The results for the Z3 gauge-Higgs model were presented in detail in our recent paper [2]. From 
a physics point of view the interesting aspect of the analysis of the dual model is the possibility 
to study a theory with gauge interactions at weak coupling for low temperature and high density. 
In such a situation one expects condensation phenomena that give rise to phase transitions as a 
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Figure 4: Lhs.: Results for (U) in the strong coupling region of the Z3 gauge-Higgs model (77 =0.1 and 
P =0.6) versus the chemical potential jl. Rhs.: The corresponding occupation numbers for temporal and 
spatial plaquettes (P t ,P s ), for spatial fluxes (S), and for positive and negative temporal fluxes (7+, 71) . 



function of the chemical potential (see [5] for examples of related condensation phenomena at 
strong coupling or with non-gauge interactions). In Fig. 4 we provide an example for such a 
condensation: The lhs. plot shows the plaquette expectation value (U) as a function of ji, which 
obviously undergoes a pronounced first order transition near ji c ~ 2.6. The rhs. plot shows the 
corresponding average plaquette and flux occupation numbers, illustrating that the transition can 
be understood as the condensation of fluxes and plaquette occupation numbers. 

Even more interesting is the physics of the U(l) gauge-Higgs model, which we currently study 




Figure 5: In the lhs. plot we show a 3-D plot of the plaquette expectation value (U) as function of = 
K*( 2 ) = k and j3 for =0 and jl^ = 0. The contour lines in the base of the lhs. plot and on the rhs. are 
equally spaced lines of constant (U). 
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Figure 6: Plaquette expectation value (U) (lhs.) and plaquette suscpetibility %u (rhs.) as function of the 
chemical potential (8 3 x 100, k = 6.4, j3 = 0.6, A = 0). 

for = A( 2 ) =0 and degenerate masses m^ 2 \ i.e., = K*. At K" = which 

corresponds to pure U(l) gauge theory, the system undergoes a first order transition at j8 c ~ 1, 
which is signalled by a pronounced disccontinuity in (U). This transition persists also for the two 
flavor model at finite K and in Fig. 5 we show our results from the dual simulation for (U) as a 
function of /3 and K at ji^ = [1^ = 0. One clearly observes the transition that separates the strong 
coupling phase from the weak coupling phase with <C/> ~ 1 

In Fig. 6 we plot the plaquette expectation value (U) and the corresponding susceptibility Xu 
as functions of the chemical potential ji for the parameter set K = 6.4, j8 = 0.6 (strong coupling 
regime). We find that the observables are independent of the chemical potential up to \i c ~ 2.9 
(Silver Blaze phenomenon), where the system undergoes a pronounced first order transition, which 
as in the Z3 case corresponds to a condensation of the dual degrees of freedom. 
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